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Quantum algorithms are conventionally formulated for implementation on a single 
system of qubits amenable to projective measurements. However, in expectation 
value quantum computation, such as nuclear magnetic resonance realizations, the 
computer consists of an ensemble of identical qubit-systems amenable only to ex- 
pectation value measurements. The prevalent strategy in such expectation value 
implementations of quantum algorithms has been to retain the conventional for- 
mulation's unitary operations but modify its initialization and measurement steps 
appropriately. This naive approach is not optimal for Grover's algorithm and a 
shortened version for expectation value quantum computers is presented. 



Quantum algorithms EJ are usually formulated with the idea that each will 
be implemented on single quantum system amenable to projective measure- 
ments. In contrast, most experimental implementations have used room tem- 
perature solution state NMR, which involves an ensemble of non-interacting 
identical quantum computers and for which the_miLy measurement outcomes 
available are expectation values of observables na n. Consequently quantum 
algorithms have been adapted for such expectation value imantum computers 
by modifying the preparation and measurement stages aU. However, these 
adaptations have been constructed so that the expectation value quantum 
computer mimics its single quantum system relative as closely as possible. As 
such, they do not use the multiple identical quantum computers in the ensem- 
ble advantageously. In this article we show how the ensemble can be used to 
shorten Grover's search algorithm for expectation value quantum computers 
(a more detailed description appears elsewhere u) . 

We consider an expectation value quantum computer, consisting of an 
ensemble of identical, non- interacting quantum jComputcrs, which (i) can be 
prepared in any pure or pseudopure initial statccrEl, (ii) to which any unitary 
transformation can be applied, and (iii) for which the expectation value (EV) 
of (Tz{k) can be measured for each qubit (labeled by k). This is the case with 
current NMR realizations. 

The standard translation of a quantum algorithm to a form suitable for 
an expectation value quantum computer can be illustrated via Grover's al- 
gorithm for searching a database containing a single marked item in one of 
N possible locations (for convenience assume that N = 2^ for some inte- 
ger L; this is not a severe restriction). The database will be represented as 
X — {0,1, . . . , N — 1} and the search proceeds with the aid of the following 
oracle, invoked at unit cost: f{x) ~ if x =^ s and f{x) = 1 if x ~ s where 
s denotes the marked item's location. A classical sequential search requires 
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N/2 oracleuiueries on average to locate the marked item correctly. A quantum 
algorithm ExLl uses an L qubit data register and an oracle unitary operation, 
defined on the computational basis states as tlf \x) :— (— l)-'^^^^ \x) and ex- 
tended linearly to superpositions of these. The data register is initialized to 
\ipi) l/\/2^X]x=o^ l'^)' followed by repeated applications of the Grover it- 
erate G := DUf. Here D (j2x=o l^^)) Ylx=o + 2 (c)) \x) , where 
(c) — c-x/N . After to applications of G, the data register's state is i 



cos [(2m + l)(?/2l v-^ , > , , 

= sin [(2to + 1)0/2] |,s) + ' ^ Y.\x) (1) 



where cos0 = 1 — 2/N. After [7r/20j applications of G the system's state 
is approximately |V') = |s); a projective measurement (PM) in the compu- 
tational basis would yield s with probability of error, which we henceforth 
ignore, at most For N ^ 1, 6 ^ 2^J\/N and the algorithm requires 

about \y/N 7r/4j oracle invocations to locate the marked item with near 
certainty, giving a quadratic speedup over the classical sequential search. 
To convert between PM and EV outcomes note that, when = |s), 
{(Jz{k)) := (V^l crz(fc) IV') = (—I)'** where s = s^.-.si is the binary repre- 
sentation of s. This reveals the marked item if it is possible to distinguish 
between {(Jz{k)) = -f-l and {az{k)) — —1 for each qubit. This standard EV 
version of Grover's algorithm has been used in all NMR realizations to date. 

What is crucial for the success of the above translation, is that the sign 
of {<Jz{k)) for each qubit can be determined. This depends only on which PM 
outcome would be more likely: positive or negative for PM outcomes or 1 
respectively. Therefore, for expectation value quantum computers, it is not 
essential that the state prior to measurement be |s), but rather one for which 
the most likely PM outcome is s. This can be attained with to < \\/N 7r/4j 
applications of G, in which case, Eq. (0) yields 

(a,(fc)) = {-ly-Am (2) 

where the EV attenuation 

8m^[{2m + l)9/2]N-l 

N-1 ^ ' 

determines the magnitude of the EV. The EV attenuation increases mono- 
tonically with with respect to to from prior to any applications of G to 
approximately 1 in the case of the standard EV version. This leads to a trun- 
cated EV version of Grover's algorithm in which is identical to the standard 
EV version except that it terminates after the minimum number of applications 
of G such that it is still possible to distinguish reliably between (^^(fc)) > 
and (o'z(fc)) < for all data register qubits. In an ideal case (i.e. an infinitely 
large ensemble and infinitely precise measurements) the truncated EV version 
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succeeds after one application of G. In practice, noise and statistical limita- 
tions effectively establish a threshold EV attenuation, Ath, below which the 
sign of the EV cannot be discerned reliably; the algorithm then succeeds when 
Am > Ath- The threshold will depend on experimental details but, for many 
practical expectation value quantum computers it will possible to truncate 
the algorithm compared to the standard EV version. 

Searching a database containing more than one marked item is compli- 
cated by the fact that the standard version of Grover's algorithm is not deter- 
ministic. Here, after about [y^N/AI 7r/4j applications of G, the data register 
is approximately in the state lipj) = I/n/^X^kss I^) where M is the number 
of marked items and S denotes the set their locations El. A PM will yield a 
location of one of the marked items with certainty. However, the bitwise EV 
technique described earlier may fail in some cases since averaging over mope 
than one possible database location can produce devastating cancellations 0. 
One way to avoid this and still incorporate the truncation described above is 
to use a filtered EV technique based on iterated runs of the algorithm. The 
iteration step produces the fc + 1 bit of a marked item's location if the first 
k bits have been determined to be Sk, Sk-i, • ■ ■ , si by combining the EVs of 
two runs of the algorithm. In the first the algorithm is run without modifica- 
tion, while in the second a correlation operation Gk+i{sk, . . . , si) is inserted 
between the last application of G and the EV measurements, where 

(7fc+i(sfc, ...,.s,)\xL... xi) (a,(fc + l))^^---^) . . . x^) , (4) 

with 

g{xk ...xi) := {xk ® Sfc ® l){xk-i ® Sk-i ® 1) . . . (xi ® Sa ® 1) 1. (5) 
Averaging the EVs obtained from the two runs gives i 

xeS' 

identical to an EV performed on ^/VM^^^g, \x) where S' :— {x £ S : Xk = 
Ski ■ ■ ■ ,xi — si}, and whose sign gives Sk+i for some marked item whose first 
k bits are Sfc, Sfc-i, . . . , si. The process is initiated by inspecting the sign of 
(cz(l)) , giving si for a marked item. The entire filtering procedure requires 
log2 N runs of the algorithm and 0((log2 N)"^) operations to compute g for all 
stages. Note that the amplitude of the EV scales as 1/M and thus for a given 
experimental realization there will be a threshold (in M) beyond which this 
scheme fails. We only consider cases where M lies beneath this threshold. 

The truncated EV version of Grover's algorithm for one marked item is 
extended to multiple marked items if the filtering scheme described above is 
used. After m applications of G, 
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where S' includes filtering conditions and the EV attenuation is 

_ sin^ [{2m+l)0/2]N~M 

N-M ' ^ ' 

and the algorithm can be terminated whenever Am > Ath- The reduction 
in the number of applications of G depends on the actual threshold EV 
attenuation, Ath, and the maximum tolerable threshold EV attenuation, 
^stand = ^/M, for the standard EV version □. The minimum number of 
applications in a truncated EV version is 

2 



TOtrunc ~ mstand - ErCSin V^thMstand + (1 " AhMstand)^^ (9) 
TT 

where rTiptand is the number of applications of G required for the standard EV 
version u. 

Whenever a quantum computing device only offers expectation value out- 
comes, then the truncated version of Grover's algorithm is superior to the 
standard version. Here the ensemble's ability to produce decisive expectation 
values may be used advantageously and a naive translation from the single 
quantum system version of an algorithm is not always optimal. The extent 
to which this is applicable to other algorithms is not clear and would be a 
worthwhile issue to investigate. 
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